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Differential cross sections for 3 He-a scattering were measured in the energy range 
up to 3 MeV. These data together with other available experimental results for 3 He 



q ' +a and 3 H +q scattering were analyzed in the framework of the optical model using 

double-folded potentials. The optical potentials obtained were used to calculate 
the astrophysical S-factors of the capture reactions 3 He(a,7) 7 Be and 3 H(a,7) 7 Li, 
and the branching ratios for the transitions into the two final 7 Be and Li bound 
states, respectively. For 3 He(a, 7) 7 Be excellent agreement between calculated and 
experimental data is obtained. For 3 H(a, 7) 7 Li a 5(0) value has been found which 
is a factor of about 1.5 larger than the adopted value. For both capture reactions a 
similar branching ratio of R = (7(71) / '(7(70) ~ 0.43 has been obtained. 
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I. INTRODUCTION 



The reaction 3 He(a,7) 7 Be determines together with other reactions the branching ratio 
between the ppl and (ppll + ppHI) chain in hydrogen burning of main-sequence stars. The 
magnitude of the cross section of this capture reaction is of special interest for the solar 
neutrino problem The mirror reaction 3 H(a,7) 7 Li is the main source for 7 Li production 
in primordial nucleosynthesis 0. 

Experimental data for the 3 He(a,7) Be cross section at subCoulomb energies have first 
been obtained by Parker and Kavanagh ||. Further experiments observing the capture 
7-rays have been performed by Nagatani et al. Q, Krawinkel et al. f5j, Osborne et al. |§, 
Alexander et al. 0, and Hilgemeier et al. [§]. Capture cross sections observing the decay 
of the 7 Be residual nucleus have been measured by Osborne et al. M, Robertson et al. 0, 



Volk et al. jKj and Hilgemeier et al. ||. Theoretically, Tombrello and Parker [|TT|] have first 
succeeded in describing the energy dependence of the astrophysical S-factor very well in a 
direct-capture model. Further calculations in the framework of the potential model have 
been carried out by Kim et al. [0 and Buck et al. |13| |. Analyses using microscopic theories 
based on the resonating group method have been performed by Walliser et al. [I4j, Kajino 
and Arima |15[], Mertelmeier and Hofmann ||16|| , Langanke [|I7[], Kajino ]18[| and Liu et al. 
P~5[ . The adopted value for the astrophysical S-factor chosen by Bahcall [|TJ to calculate the 
expected solar neutrino flux is S(0) = (0.54 ± 0.03) keV b. 

In the case of the 3 H(a,7) 7 Li reaction three sets of experimental data at subCoulomb 
energies are published by Griffith et al. [20|1, Schroder et al. ||21|| , and Burzynski et al. 
Several microscopic calculations have been performed in order to analyze this reaction 



||l^ , p!6|Jl7| , p!qj2^ , |24| . The presently adopted S-factor value used in standard hot big bang 
model studies is given by S(0) = 0.064 keV b |25 . 



In this work we analyze the experimental data of both capture reactions, 3 He(a,7) 7 Be 
and 3 H(a,7) 7 Li, for energies Eqm < 1-4 MeV and 0.6 MeV, respectively. The calculations 
were performed in the framework of the direct-capture model. The most important ingre- 
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dients in this model are the optical potentials for the bound and scattering states. These 
potentials are determined using the folding procedure. The strengths of these potentials are 
adjusted to the experimental scattering data. Therefore, we measured the differential cross 
sections for 3 He - 4 He elastic scattering in the range ( 3 He) < 3 MeV. 

In the next section we describe the direct capture model and the folding procedure for 



the optical and bound state potentials. In Sec. |TTT| we present the 3 He - 4 He and 3 H - 
4 He cluster potentials derived from elastic scattering measurements. Finally, in Sec. [IV| the 
results for the astrophysical S-factors of the capture reactions are given and compared with 
the experimental data. A summary is given in Sec. |V|. 

II. DIRECT CAPTURE MODEL AND FOLDING PROCEDURE 

Potential models are based on the description of the dynamics of nuclear processes by a 
Schrodinger equation with local potentials in the entrance and exit channels. Such models 
are the Optical Model (OM) for elastic scattering, the Distorted- Wave Born Approximation 
(DWBA) for transfer and the Direct Capture Model (DC) for direct capture reactions. 

The most important ingredients in the potential models are the wave functions for the 
scattering and bound states in the entrance and exit channels. In calculations performed by 
our group the potentials are determined by using the folding procedure. In this approach 
the number of open parameters is reduced considerably compared to more phenomenological 
potentials (e.g. Saxon- Woods potentials). The nuclear densities are derived from nuclear 



charge distributions |26| and folded with an energy and density dependent nucleon-nucleon 



(NN) interaction v cS 0,|8],|29 



V(R) = XV F (R) = X j J p a (ri)p A (r2)v c g(E } p a ,p A ,s)dfidf*2 ■ (1) 
The variable s in the NN interaction term is given by 

s = \R + f 2 -fi\ (2) 
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with R being the separation of the centers of mass of the two colliding nuclei. The nor- 
malization factor A accounts for Pauli repulsion effects and dispersive parts in the potential 
V(R) which are not included in the folding potential Vp(R). This parameter can be ad- 
justed to elastic scattering data and/or to bound and resonant state energies of nuclear 
cluster states. At the low energies considered in the nucleosynthesis often the imaginary 
term in the potential can be neglected. Therefore, in the potential model combined with 
the folding procedure for the potential the reaction cross sections can be calculated in many 
cases without any free parameter. 

The DC cross section is given by |Kj 



dn 

dn 2 (h) fe) 2ihi 2s^~i m J; Mb , 1 T ™'"-° |2 ' (3) 

The quantities I a, Ib and S a (Ma, M b and M a ) are the spins (magnetic quantum numbers) 
of the target nucleus A, residual nucleus B and projectile a, respectively. The reduced mass 
in the entrance channel is given by /i. The polarisation a of the electromagnetic radiation 
can be ±1. The wave number in the entrance channel and for the emitted radiation is given 
by k a and fc 7 , respectively. 

The multipole expansion of the transition matrices TM A M a M B ,a including electric dipole 
(El) and quadrupole (E2) transitions as well as magnetic dipole (Ml) transitions is given 
by 

TM A M a M B ,a = TM A M a M B ,a d& a (6) + TM A M a M B ,a ^L(^) + ^M^MaMs^ ^L(^) ■ (4) 

— * — * 

The rotation matrices depend on the angle between k a and k 1 which is denoted by 9, where 
5 = M A + M a - M B . 
Defining 



Z„_ 
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m a 


m A 
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12 \ mi m. 



we can write for the transition matrices for the electric dipole (EC = El) or quadrupole (EC 
= E2) transition 

Tlf AMa M B ,a = E %la & S a M a \ j a M a ) (j b M B - M A I A M A I I B M B ) 
x (C5j b M B -M A \j a M a )C(EC)lJ b j b 

x (l b 0C0 | l a 0) W(Cl b j a S a ; l a j b ) lF b f blB;la3a . (7) 

In the above expressions Z a , Z A and m a , m A are the charge and mass numbers of the 
projectile a and target nucleus A, respectively. The quantum numbers for the channel spin 
in the entrance channel and for the transferred angular momentum are denoted by j a and 
j b , respectively. 

For magnetic dipole transitions (MX = Ml) we obtain 

TM C A M a M B ,a = E *'° ° { & 0S « M « Ja M a ) (j b M B - M A I A M A \I B M B ) 

laja y 



x (1 5j b M B —M A | j a M a ) 

( 



X 



H\ % + ) i b J b Jl a (!a + l) W(l U« S a] l a j b ) 

m\ m\) v 



+ 2n a (-l)*>->- S ajb ^/S a (S a + 1) W(l S a3a l a , S a j b ) 
- (l a 0S a M a \ Ja M a )( 3a M a I A M B -M a | I B M B ) 
x (I A M B — M a 1 5 | J A M A ) 



x 



^ S jajb y/(I A + I)/ 1 A \ \ I <U Z a /J5^ ;Wa , (8) 



2m p c 



where W is the Racah coefficient, the are the magnetic moments and m p is the mass of 
the proton. 

The overlap integrals in Eqs. (7) and (8) are given as 

'k&isM. =J dr WO ° EC ( r ) Xia^r) (9) 
for the electric dipole (EC = El) or quadrupole (EC = E2) transition, and by 



rMl 

hjblB \laja 



dru NLJ (r) M \r) XUa (r) 



(10) 



for the magnetic dipole transition (MX = Ml). 

The radial part of the bound state wave function in the exit channel and the scattering 
wave function in the entrance channel is given by unlj{^) and Xi a ja{ r )-> respectively. The 
radial parts of the electromagnetic multipole operators are 

1 



M1 (r) 



2p 



P 3 
15 



[sin p + p cos p] , 
(p 2 — 2) sin p + 2p cos p 



<D^(r) = — [{hp 1 - 12) sinp + (12 - p 2 )p cos p 



(11) 
(12) 

(13) 



In the long wave-length approximation - applicable in our case, since p = k^r C 1- these 
quantities reduce to 



O m (r)~l , 
E1 ( r )~r , 
C E2 (r) ~ r 2 . 



(14) 
(15) 
(16) 



III. 3 He 4 He AND 3 H 4 He CLUSTER POTENTIALS 

At incident energies _Ei a b( 3 He) ^3 — 10 MeV differential cross sections for the elastic 3 He 
4 He scattering are well known |32] , |33|j3^ , |35|j36]| . In order to obtain cross section data at very 



low energies we have measured these scattering processes at 20 energies in the range from 
1 to 3.3 MeV using the windowless, differential pumped and recirculating gas target system 
RHINOCEROS installed at the Stuttgart Dynamitron accelerator J37[. In our experiment 



we used the jet configuration in which a supersonic jet produced by a laval nozzle serves as 
a nearly pointlike target zone with high density ]3"8|] . The small and fixed size of this zone 
allows a good determination of angular distributions. 
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The target zone is centered in a 50 cm-diameter scattering chamber. The detector system 
consisted of ten surface-barrier detectors mounted at fixed positions. 3 He-beam intensities 
between 5 and 30 /zA were used. In order to normalize the data (assuming that the 3 He- 20 Ne 
scattering is dominated in this energy range by Rutherford interaction), a small quantity of 
20 Ne gas was admixed to the 4 He gas in the jet. 

The experimental differential cross sections for the 3 He - 4 He scattering are shown to- 
gether with older data of Barnard et al. ||3^] and Chuang |34| in Fig. |l|. The results of different 



phase-shift analyses H33| , ^5l , |36| for projectile energies _Ei a b( 3 He) > 3 MeV are presented in 
Fig. |. 

We have calculated the differential cross sections and phase shifts in the framework of 
the OM. For the calculation of the real part of the optical 3 He - 4 He potential we used the 
folding procedure as described in Sec. ||. The folding potentials (Eq. 1) were determined 
using the computer code DFOLD [Q. The imaginary part was neglected because the flux 
into other channels is very small. Together with the spin-orbit term, the optical potential is 
given by 

V(R) = \V F (R) + A so i^^L • s (17) 

with a spin-orbit normalization factor A so . 

As result of a fit to the experimental data given in Figs. p] and [| we obtain a parity- 
dependent potential. The normalization factors A together with the volume integrals per 
nucleon pair Jr are listed in the upper part of Table |. The spin-orbit normalization factor 
A so = —0.162 fm 2 has been determined from the splitting of the phase shifts for the L = 3 
doublet. The agreement between the experimental data and the results of the OM calculation 
is excellent in the whole energy range up to 14 MeV as can be seen in Figs. [I] and 0. 

For the 3 H - 4 He system phase-shift analyses of experimental cross section data have been 
performed in the energy range _Ei a b( 3 H) =3 — 10 MeV J33|. We have calculated these phase- 



shifts in the OM. The optical potential is again determined using the folding procedure. The 
result of our OM fit to the phase-shift data is shown in Fig. 0. The normalization factors 



A together with the volume integrals per nucleoli pair Jr are presented in the lower part of 
Table |. For the spin-orbit normalization factor A so now we obtain A so = —0.136 fm 2 . Again 
the agreement between the experimental and calculated data is satisfactory (s. Fig. |3|). 

The volume integrals of the parity-dependent potentials for a-scattering on the mirror 
nuclei 3 He and 3 H only differ by about 0.5% and 3% for the even and odd partial waves, 
respectively. The values of the volume integrals for the even partial waves in 3 He - a and 3 H 

- a scattering is comparable with the value Jr = 445.7 MeV fm 3 obtained in the analysis 
of a - a scattering using the folding procedure [4(J. Furthermore, the values of Jr for 3 He 

- a scattering are compatible with the results of a systematic analysis of 3 He-scattering on 
several nuclei ||41|J42||. 



In a next step we used the double-folded potential as a suitable cluster-cluster poten- 
tial and calculated bound states and single-particle (single-cluster) resonances. The wave 
function Ujvxj(r) which describes the relative motion of the respective 3 He - a and 3 H - a 
system is characterized by the node number N and the orbital angular momentum L. The 
N and L values are related to the corresponding quantum numbers and /j of the three 
nucleons forming the 3 He and 3 H cluster, respectively: 

3 3 

Q = 2N + L = Y / 2n i + h = ]>> = 3 . (18) 

i=l i=l 

Thus for both systems 7 Be = 3 He ® a and 7 Li = 3 H <g> a one expects two cluster states 
with L = 1 (N = 1) and L = 3 (N = 0). Both states split into doublets with J = L ± 1/2 
because of the spin-orbit potential resulting from the motion of the A = 3 particle with spin 
1/2 in the field of the a-particle. 

In the calculations of the cluster states the centroid energies of both the bound {L = 
1 : J n = 3/2-,l/2~) and quasi-bound (L = 3 : J w = 7/2", 5/2") states for 7 Be and 7 Li 
(s. Fig. |5|) are reproduced by the central parts of the 3 He - a and 3 H - a odd potential, 
respectively. The splitting of the energies of the quasi-bound state doublets (L = 3) in 
7 Be and 7 Li is reproduced by the spin-orbit potential as determined by the OM calculation. 
However, for the bound-state doublets (L = 1) in 7 Be and 7 Li smaller spin-orbit potentials 
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are necessary. For the energy splitting of the bound-state doublets a spin-orbit normalization 
factor A so = —0.07 fm 2 has to be used for both nuclei 7 Be and 7 Li. 

As a further test of our folding potential we calculated the charge distribution of the 7 Li 
nucleus (in its ground state) by folding the experimental charge distribution of 3 H and 4 He, 
which we have already used in our double-folding procedure, with the radial wave function 
wjv=i,i=i,j=3/2 and by assuming a spherical shape for the folded distribution. The result of 
this calculation is shown in Fig. |] together with the experimental charge distribution, as 
measured by electron scattering on 7 Li. It can be seen that the calculation overestimates the 
experimental density in the nuclear interior, whereas in the surface region, near r = 2.5 fm, 
the experimental values are slightly underestimated. But the rms radii of both distributions 
are almost identical: < r 2 > 1//2 ^; 2.40 fm. 



IV. CAPTURE REACTIONS 

In Fig. |5] a schematic presentation of the direct capture processes 3 He(a,7) 7 Be and 
3 H(a;,7) 7 Li is given. In the low energy range capture transitions can only occur into the 
ground and first excited states of 7 Be and 7 Li, respectively. 

The theoretical cross section cr th is obtained from the DC cross section <t dc , given in Eq. 
3, as sum over both final states i = 1, 2 by 

^ = £<^ D ° . (19) 

i 

The computation of the cross section a th was performed using the computer code 
TEDCA ||43| . As input three data sets are necessary: (i) isospin Clebsch-Gordan coeffi- 
cients, (ii) spectroscopic factors Si which specify the cluster probability of the final states 
7 Be = 3 He (g> a and 7 Li = 3 H eg) a, respectively, and (iii) optical potentials for the calcula- 
tion of the wave functions in the entrance and exit channel. In our case the Clebsch-Gordan 
coefficients are C\ = 1. The spectroscopic factors Sj have been taken from the work of 



Kurath and Millener [44 . The numerical values are given as Si = 1.174 and S2 = 1.175. 



In order to calculate the bound state wave function in the exit channel and the scattering 
wave function in the entrance channel (Eqs. 9,10), the folded potentials are used which have 
already been determined in Sees. [H] and [[II]. That means that all the necessary information 
for the calculation of the DC reaction is known and no parameter has to be adjusted to the 
experimental capture reaction data. 

In Fig. P the experimental values of the astrophysical 5-factor for the reactions 
3 He(a,7) 7 Be |,||| and 3 H(a,7) 7 Li |20|,|2ll are shown together with the results of DC cal- 
culations using parity-dependent folding potentials. The experimental data of Krawinkel 
et al. || are renormalized by a factor of 1.4 as suggested by Hilgemeier et al. |§. 

For the 3 He(a,7) 7 Be reaction the agreement between the experimental and calculated 
data is excellent. A linear extrapolation for E — > gives S'(O) = 0.516 keV b and S(0) = 
—3.67 • 10~ 4 b (E in keV). The values of S(0) agree excellently with the experimentally 
determined S(0) factor, for which || gives a weighted average of S(0) = (0.51 ± 0.02) keV 
b, and with the adopted value |1[ of S(0) = (0.54 ± 0.03) keV b. The calculated branching 
ratio R for the transitions to the first excited state and the ground state has the value 
R = 0.43. This value was found to be nearly energy-independent in the energy range - 1.4 
MeV and agrees well with both the experimental data |§|3|U and the results of microscopic 
calculations |24| . 



The two presently available measurements of the low-energy 3 H(a,7) 7 Li reaction differ 
from each other by roughly 30% in total magnitude as well as in their determination of the 
branching ratios for the transitions into the two final bound states. The calculation within 
the potential model gives results for the absolute magnitude of the S-factor which favour 
rather the older experimental data of Griffith et al. p0| . The calculated branching ratio 



is again nearly energy-independent and has likewise a value of R = 0.43. This branching 
ratio, however, is noticeably larger than the experimentally observed value R = 0.32 ± 0.01 



21j| . Our value agrees with the data of [EDJ which give an energy-independent average of 



R = 0.43 [|24}|. The above considerations are comparable with the results of microscopic 



RGM calculations [ 24| . A linear extrapolation for E — > gives S(0) = 0.100 keV b and 

10 



5(0) = — 1.02- 10~ 4 b (E in keV). The values of 5(0) are nearly twice as large as the adopted 
value of 0.064 keV b ||25|| , obtained from an energy-independent extrapolation, and somewhat 
smaller than the extrapolated value 0.14 keV b extracted from the data of Schroder et al. 



21]. However, it is in excellent agreement with the results of different RGM calculations 



1J,[T3.|TJ.|23|J24| giving a value of 5(0) « 0.1 keV b. 

In Fig. [7] the multipole contributions for both capture reactions are shown. The main 
contribution is the El transition. Because of the missing centrifugal barrier, the DC transi- 
tions from the s-wave in the entrance channel to the final L = 1 states are dominating. For 
the 3 H(a,7) 7 Li reaction the contributions of the higher partial waves for the total S-factor 
can even be neglected. As can be seen on the right-hand side of Fig. [5] the curve of the total 
S-factor is almost identical to the s-wave contribution. The Coulomb barrier is lower for 
3 H(«,7) 7 Li than for the mirror reaction. Therefore, the influence of the centrifugal barrier 
is more pronounced. 

As already discussed in Sec. [HI], a parity-dependent optical potential is necessary to 
describe the experimental scattering data. We calculated the astrophysical S-factor with a 
parity-independent potential using the A-parameter for the dominating even partial waves 
given in Table |I[ As shown in Fig. [|, it is impossible to reproduce the experimental data 
with such a parity- independent potential. The enhancement for energies E > 0.7 MeV is 
due to the Ml-contribution of the p-wave. 

We also tested the sensitivity of the astrophysical S-factor by changing the strength of the 
optical potentials in the entrance channel. A variation of the potential depth of ±1% leads to 
an energy-independent change of the S-factor of ±2%. The S-factor increases with growing 
potential depth. The reason for this behaviour is that with growing nuclear potential depth 
the Coulomb barrier becomes smaller and therefore the penetration probability is enhanced. 
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V. SUMMARY 



Differential cross sections for the elastic scattering of 3 He particles on 4 He and phase 
shifts of both 3 He - 4 He and 3 H - 4 He scattering have been analyzed up to energies of 
about 10 MeV in the framework of the optical model. The potential was deduced by a 
double-folding procedure using a density-dependent effective nucleon-nucleon interaction. 
The experimental data are described satisfactorily by this optical-model calculation. 

Using the double-folded 3 He - a and 3 H - a potentials as cluster-cluster potentials, 
we calculate the bound and quasi-bound doublet states in 7 Be and 7 Li, respectively. The 
excitation energies of these states as well as the charge distribution of 7 Li are well reproduced 
in this potential model. 

The optical potentials obtained from the fit to the elastic scattering data have been used 
to calculate the astrophysical S-factors of 3 He(a,7) 7 Be and 3 H(a,7) 7 Li within the direct 
capture model. Using this method no parameter has to be adjusted to the experimental 
reaction data. 

In the case of the reaction 3 He(a,7) 7 Be we obtain S(0) = 0.516 keV b. This value is 
in excellent agreement with both the average of the experimental data || and the adopted 
value |IJ. The branching ratio for the transitions to the first excited state and the ground 
state results in R = 0.43, likewise in very good agreement with the known experimental and 
theoretical data. 

For the reaction 3 H(a,7) 7 Li the three presently available measurements differ from each 
other by roughly 30% in total magnitude as well as in their determination of the branching 
ratio for the transitions into the two final 7 Li bound states. In agreement with calculations 
which have been done in the framework of the Resonating Group method or on the basis of a 
microscopic potential model, our calculations predict the S'-factor to increase with decreasing 
energy resulting in S(0) = 0.10 keV b. This value is a factor of about 1.5 larger than the 
adopted value p5| . For the branching ratio we obtain R = 0.43. This value is in good 
agreement with some theoretical results, but is not compatible with the recently measured 



12 



value of R = 0.32 ± 0.01 
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TABLES 

TABLE I. Normalization factors A and volume integrals per nucleon Jr of the optical potentials 

partial wave A jR(MeV/fm 3 ) 

3 He - a even (s,d) 1.452 469.0 

odd (p,f) 1.844 595.6 

3 H - a even (s,d) 1.525 466.8 

odd (p,f) 1.890 578.5 
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FIGURES 

FIG. 1. Experimental differential cross section for the elastic 3 He - 4 He scattering for projectile 
energies between £'i a b( 3 He) = 1.2 and 3.0 MeV together with data of Barnard et al. |}2| and Chuang 
The solid lines are the result of the present OM calculation, the dashed lines give Rutherford 
scattering. 

FIG. 2. Phase shifts deduced from experimental 3 He - 4 He scattering data given by Spiger and 
Tombrello (A, @), Boykin et al. (o, §§) and Hardy et al. (n, f3§). The solid lines are the result 



of the present OM calculation. 

FIG. 3. Phase shifts deduced from experimental 3 H - 4 He scattering data given by Spiger and 
Tombrello [|33||. The solid lines are the result of the present OM calculations. 



FIG. 4. Comparison of the experimental charge distribution of 7 Li (2(| (dashed line) with the 
distribution calculated in the potential model (solid line). 

FIG. 5. Schematic presentation of the level scheme for the reactions 3 He(a,7) 7 Be and 
3 H(a, 7 ) 7 Li. 

FIG. 6. Calculated astrophysical S-factor using the potential model compared with the exper- 
imental data for 3 He(a,7) Be (upper part: closed circles ||, open circles |J, triangles ||), and 



H(a, 7) Li (lower part: open circles [£Cj], closed circles [21], triangles |22||). 



FIG. 7. Multipole contributions to the astrophysical S-factors for the reaction 3 He(a,7) 7 Be 
(left-hand side) and 3 H(a,7) 7 Li (right-hand side). The experimental data are the same as in 
Fig.§. 

FIG. 8. S-factor for 3 He(a, 7) 7 Be calculated with a parity-independent potential in the entrance 
channel. The experimental data are the same as in Fig. |6|. 
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